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Perturbation  Theory  for  Implosions 

Cathleen  S,  Morawetz 

1.   Introduction. 

If  a  piston  is  pushed  into  a  closed  tube  of  gas  a  shock 
is  created  which  races  ahead  of  the  piston  to  the  end  of  the 
tube.   There  it  is  reflected  and  returns  back  to  the  piston. 
This  process  is  very  well-known  and  the  flow  has  been  com- 
pletely determined  for  arbitrary  piston  motions.   In  this 
paper  ^^^e  are  concerned  with  an  analogous  problem  in  three 
dimensions,  the  so-called  implosion  phenomenon,   A  spherical 
shell,  or  "piston",  is  contracted  toward  its  center,  thereby 
creating  a  shock  which  precedes  it  into  the  center.   The  prob- 
lem is  to  determine  the  flow  inside  the  shell  and  to  decide 
whether  the  shock  is  reflected  out  from  the  center, 

Guderley  [1]  has  determined  such  a  flow  for  a  special 
choice  of  piston  velocity.   In  this  particular  case,  the  shock 
is  reflected  at  the  center  and  expands  bad':  toward  the  shell 
jUst  as  in  the  one-dimensional  case.   One  very  interesting 
result  is  that  the  pressure  behind  the  reflected  wave  is  about 
twenty-six  times  the  incident  pressure  while  in  the  one- 
dimensional  case  the  ratio  is  only  eight,   VJe  shall  call  this 
flow  the  basic  flow. 

The  question  arises  immediately  as  to  whether  this  flow 
pattern  is  general  or  just  a  special  effect,  resulting  from 
the  choice  of  piston  velocity,  that  happens  to  resemble  the 
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one-dimensional  case.  We  would  like  in  fact  to  show  that  the 
flow  resulting  from  an  arbitrary  piston  motion  behaves  like 
the  basic  flow  near  the  center  and  that  therefore  the  reflec- 
tion pattern  at  the  center  is  essentially  the  same  in  all  cases. 

By  using  the  method  of  perturbation  it  is  possible  to 
show  that  if  the  piston  is  moved  with  a  velocity  somex>fhat 
different  from  that  which  yields  the  basic  flow  the  differ- 
ence between  the  two  flows  disappears  at  the  center.   Thus 
it  is  to  be  expected  that  pistons  contracting  v;ith  arbitrary 
velocity  produce  flows  with  reflection  properties  similar  to 
the  basic  flow. 

In  addition  It  is  shovm  how  to  determine  the  perturb- 
ation terms  that  are  due  to  a  shift  in  the  time  scale  or  in 
the  distance  scale,  see  sections  13  and  ll;  of  Part  III, 

The  basic  flow  is  the  particular  flow  inside  a  contract- 
ing spherical  piston  which  can  be  described  by 

u  =  A"^  r  t"^  UQ(r~^t) 

c  =  \-^  r  t'^  Cq^^**^*^ 
p  =  X"^   r^  t"2  P^d-'^t) 

Z  —  A.    V     o        6 

Here   u,  c,  p   and   z  are  the  velocity,  soiond  speed, 
pressure  and  shock  velocity,  t  is  time  and  r   is  the  dist- 
ance from  the  center.   U  ,  C^  and  P^  are  all  fimctions  of 
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one  variable,  >5  =  r~ 't,  which  are  chosen  so  that  the  differ- 
ential equations  of  motion  and  the  shock  conditions  are  satis- 
fied, \   -   I.39I1.7   1"^  a  fixed  constant  and  cannot  be  chosen 
arbitrarily  as  vje  shall  see  later. 

We  shall  consider  flows  produced  by  an  appropriate  piston 
motion  which  can  be  described  by 

u  =  X"^  r  t"^  (U^(^])  +  s  U^(>'^,r)) 
c  =  X"^  r  t"^  C^(Y])(1  +  e  C-^(n,r)) 

p  =  X"2  r^  t"^  Pq^*')^^^  ^  ^  Pl(ij,r)) 

-X  2 

vjith  M  =  r"  t.   Here   e   Is  small  and  terms  of  the  order   e 

are  neglected.   In  this  case  a  preliminary  but  lengthy  calcul- 
ation indicates  that  if  vie   keep  n  constant  and  let   r  — i>  0 

^^/e  find  that  the  quantities   U,  ,  C,   and  P,   all  tend  to  zero, 

2 
approximately  like   r  , 

At  some  future  date,  it  is  Intended  that  these  calcul- 
ations will  be  extended  in  a  more  rigorous  way. 

Thus  the  flow  behaves  like  the  basic  flow  near  the  center. 
Furthermore,  the  singularities  that  occur  at  the  center  of  the 
flow  vjlll  be  the  same  as  those  at  the  center  of  the  basic  flow. 
Finally,  the  shock  will  be  reflected  at  the  center  and  large 
pressures,  similar  to  those  in  the  basic   flow,  ^^flll  occur. 
Part  I ,   Diff e rential  Equations  and  Boundary  Conditions , 
2  Formulation, 

V/e  assume  that  the  spherical  piston  Is  contracting  with 
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velocity  u   on  the  enclosed  gas  which  is  originally  in  a 
constant  state  at  rest,  that  is,  it  has  velocity  u   =  0, 
pressure  p   =  constant   and  densa.ty   p   =  constant.   The 
velocity  u,  pressure   p,  and  density   o   resulting  in  the 
gas  satisfy  the  equations  for  spherically  syrametric  motion, 
(see  [2],  p.  ]+l8), 

u,  +UU   +-P   =0 
t      r   p  ^r 

(1)  p,  +  u  p   +  d(u   +  — )  =  0 
^              jt     jT       j   r    r 

(pD~^)^  +  u(pp"^)^  =  Oc 

It  is  assumed  here  that  the  gas  is  polytropic. 

The  given  path  of  the  piston  is  desci'ibed  by 

(2)  r  =  Rp(t) 
and  at  the  shock,  the  conditions, 

p{u   -  z)  =  -p^   z 

(3)  p(u  -  Z)2  +  p  =  jD^  z^  +  p^ 

1_L^(U  -  Z)^  +  (1  -  1X^)0^  =  [l^     Z^  +  (1  -  ti,2)c^ 

hold  where   z   is  the  shock  vel&clty,  c  =  ]/ -^  is  sound 

?    Y  -  1 
speed  and  u,  ^  =  -^ — ; — ^  » 

■^     Y  +  1 

Guderley  has  found  that  if  appropriate  neu   variables  are 
introduced  special  solutions  of  (1)  can  be  found  by  solving 
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one  non-linear  differential  equation  of  first  order.   One  such 
solution  represents  a  flow  between  a  contracting  shock  and  a 
reflected  expanding  shoclc.   We  use  the  same  change  of  variables* 
We  introduce  instead  of   t   the  independent  variable 


and  the  dependent  variables 


U  =  Xtr"-"-  u 


(5)  C  =  Xtr"^  c 


P  =  X'tV^  p 


where  X  is  an  appropriate  positive  constant. 
From  (I4.)  X'je  obtain 

ir  \  dr,  _  ^  dr  ^  dt 

(6)  Y  -  -^  —  +  —  • 

Using  ()|.),  (5)  and  (6)  the  equations  of  motion  become, 
with  a  =  X"  , 

V'.p  pP 

(7)   ^U^d  -  U)  +  aUrU^  -  r'''^c/-  -^   +   ay'-^C^  -^  +  U(aU-l) 

+  Zar'-^C^   =  0 

-^jU   +  arU^  +  -^^    (1  -  u)  -  2  -^  (1  -  U)  +  aU  -^ 

rC 

-  2aU  -—■  +   3aU  =  0 
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P^  (1  -  u)  -  1^  (1  -  u)  V  ^  ^^T 


'1  -       r        2a^  x.  r 

"   Y-1  C 

Y-1 
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where    £)_=  yP/C    . 

These    equations  may  be   rewritten  In   the   form, 

rP 
Dr|U     =  A   -    a(C^   +  U(l-U))rU     -   ay'^C^  -^ 


1 


(8) 


D  ^T-*-  =  B   -    a 


111  rU     - 

2     ^  r        1-U     C 


aUD  "'"r        a(Y-l)        C^     ^^r 


2y      TiMiT    P 


Tjp  rP 

D  ^  =  E   -    ayrU     -   a(C^   +  U(l-TJ))   -^ 


Inhere 


A  =    (1   -   U)(l   -    aU)U  -    (3aTJ  -   2y~-^(1   -    a)  )C^ 


(9) 


B   =    (1   -    U)(l   -    aU)    -    (y  -1)[(1   -    alJ)    -   |    (1   -   a)]U 

-    [a   +  y"-^(1   -   a)(l   -   U)-^]c2 


E   =    (I   -   U)[2(l   +  y)(1   -    aU)    -   3y]    +  y(1   -    aU)    -    2aC' 


D  =    (1   -   IJ)2   _   c^. 


For   further  reference  we   note   here   the   identities 


-A   +  -^  B(l   -U)    +  D[3aU ^    (1   -    aU)]    =   0 


Y-1 


Y-1 


(10) 


-A   +   y"'''E(1   -   U)    +  D[3aU  -  ^   (1   -    aU)]    =   0. 
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At  the  piston  we  have 


(11)  U  =  U  =  Xv^t"'^   ~  R 

p     /     dt  p 


and  Immediately  behind  the  shock 

n(U  -  Z)  =  -jD^Z 
(12)     iJ(U  -  Z)2  +  P  =  p^Z^  +  X^rj^   ^^^"^Pq 

[x^(U  -  Z)^  +  (1  -  ti^)c2  =  l^^Z^  +  (1  -  [i^)7ryrr--   ~c 


2  2^2X.2^  2 


where  11=  yV/C^,    see  (7),  and 


(13)  Z  =  Xtr'-'-z. 


3 •   The  Basic  Solution • 

Guderley's  solutions  are  those  solutions  of  (8)  and  (12) 
which  are  independent  of  r.   In  such  a  solution  the  position 
of  the  shock  is  constant,  say   >T  =  r~  t  =  H   and  the  shock 
velocity  is  therefore  given  by   z  =  ar/t.   From  (13)  it 
follows  that   Z  =  1,   This  is  compatible  with  (1?)  only  if  we 
have  p   =  c   =  0  or.  In  other  viords,  if  the  shock  is 
infinitely  strong. 

From  (8),  v;e  find  for  this  case 


DrjU  =  A 


Uh)  D^jC^  =  BC 

f7 


DnP,,  =  EP 
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and  from  (12), 

U  =  1  -  ^x^ 


(15)  C  =  +^lI/i  +  ix^ 


P  =  (1  -  [I    }f>^* 


2 


Since  the  piston  is  contracting,  the  gas  moves  towards 
the  center,  i.e.,  u  is  negative  at  the  junction  of  the  piston 
and  the  shock.   From  (i+),  (5)  and  (111)  we  see  that   t  <  0, 
►■j  <  0,  corresponds  to  u  <  0,   Since   c  =  Y^2   >  0  we  have 
from  (5)  that   t  <  0  corresponds  to   C  <  0,   For  a  contract- 
ing piston  we  must  therefore  choose  the  negative  sign  in  the 
second  equation  of  (ll|-). 

In  order  that  equations  (II4.)  have  a  solution  satisfying 
the  initial  conditions  (15)>  that  does  not  become  singular  for 
D  =  0,  the  value  of   a  must  be  specially  chosen.   The  correct 
value  for   a   is  approximately   .717  with  y  =  I.I4-.   In  this 
case  then  A  =  B  =  E  =  0  when  D  =  0  but  the  ratios   A/D, 
B/D,  and  E/D  remain  finite. 

To  the  incident  flow  satisfying  (l!}.)  and  (15)  with 
a  =  ,717   a  reflected  shock  can  be  fitted.   The  flow  behind 
this  reflected  shock  is  easily  determined.   In.  the  incident 
flow  U  lies  between   (1  -  [j,  )   and  some  negative  value 


behind  the  reflected  shock  while   C   lies  between  -(-i/l  +  [i 
and  some  positive  value. 


j;;  J  -3V-0 J    EflV'i" 
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The  picture  of  the  flow  is  then  as  follows:   Up  to  some 

time,  t    say,  the  piston  is  fixed  at   r  =  R   where 

R^*"  t   =  H   and  the  state  in  front  of  it  is  constant.   Then 
o    o    o 

for  t  >  t   the  position  of  the  piston,  using  (5)  and  (6),  is 

given  by 

('  aU    , 
(16)      r  =  P.  ()i)  =  I  exp  \     y-2r-   ^ 

o        ' 
where   U^(',)   is  the  solution  of  (lli.)  and  (l5)  with  a  =  ,717. 
The  state  in  front  of  the  piston  is  still  constant  for   >,  <  H  , 
For  I;  >  H   and  up  to   i,  =  H    vrhich  represents  the  position 
of  the  reflected  shock  the  state  is  given  by  the  solution 
U^(*"|)>  C^(n),  Pq('|)   of  (111)  and  (15).   The  particle  paths  in 


this  region  are  given  by  the  equation 

' '  ^^   d> 
(17)      log  r  =  I  -,  ^jj     -^  +   constant, 

J   ""  o  '/ 
The  characteristics  in  the  (r,t)-plane  are  given  by 

(18)  if  =  (^^  i  °)^ 


see  [2],  p,  I4.6,  and  are  therefore  given  in  the  (>),r)-plane  by 
(19) 


'I' 
dr  _     1 


rd  ,  "  TT^iTTic;  • 


[(.,   Perturbation. 

The  present  treatment  consists  in  perturbing  about 
Guderley's  solution  and  studying  the  behavior  of  the  first 
order  terms-. 

We  assume  once  more  that  the  state  on  which  the  piston 
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contracts  has  zero  pressure", 

(20)  p   =  c^  =  0. 

o     o 

Up  to  time   t  =  t   +  et-,  ,  suppose  the  piston  is  fixed  at  the 

position  r  =  R   +  eR.,   and  at   t  =  t"  +  et",  ,  it  starts  moving 

o     1  o     1'  ° 

on  a  path  given  by 

(21)  r  =  Rp(.-j)  =  R^(r|)  +  e  R^{rj) 
where 


t ) 


(22)       R^(,,J  =   R^  exp  ( 


aU 


o  d'-j 

„      1-U  n 

H      o  ( 
o 


(Guderley  path] 


We  are  evidently  free  to  choose   t-,   and   R,  =  R,  (PI  )   since 

1        1     1   o 

any  change  in  these  quantities  can  be  absorbed  into   R    and 

^o* 

It  is  reasonable  to  expect  then  that  all  the  flow  quan- 
tities can  be  expanded  in  powers  of   e,  that  is, 

U  =  U^(i|)  +  e  U^(-|,r)  ... 

(23)  C  =  C^(.j)(l  +  eC^(r  r)  ...) 

P  =  lo(^j)(l  +  eP-^(yr)  ...) 

where   U^(n),'  ^■■,^'l^>    ^o^'P   ^^  Guderley 's  solution  of  (I'l.) 
and  (15)  in  Section  3.   'ibstituting  (23)  in  (8)  and  taking 
first  order  terms  in   e   v;e  find: 


It  can  be  shown  that  if  p    is  considered'  small  the 
resulting  disturbance  dies  out  at  the  center  like   r     ,  cf,[3]. 


':  !     ?.  j-Uj  7  ^     'j  C 
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.  1=1   Q  D"^arP 
2y   1-U  "^0  ^^  Ir 
'      o 

nP,   -  (^)°  U-,  +  (|)°  C  C,  -  YD"-^arU, 
I  In    ^D  L   1    ^D  C   o  1    '  o    Ir 

-  D"^(C^  +  U  (1-U  ) )arP. 
o  ^  o     o^    o     Ir, 


TT         /AxO  _  i  d  /An 

^^^^  ^D^u  -|  au  ^D^ 


7  =  U^(v,),  C  =  C^(K|) 


From  equations  (10)  we  see  that 


,    etc . 


^Y-1      DC         ^D^C    1-U 


(25) 


-1    ,E.         ,A, 
Y         (-)o  -    (-] 


d'G~    ^ D 'C   (1-U) 


0, 


Using    (25)    with    (2[i.),    we    obtain 


^     (nc.    )   -  ^V-(  lis  J 


Y-1 


1-U 


\ri  ^  TT  n  __TT        V  n  /"''-' -i 


>y         i^^Y-l      D'U        l-U^    'D'U  /     1 


arU 


Ir 


1-U, 


U  ^ 

o        2  „ 

1-J      Y-1  Ir 


(26)       v-^.|F,  )  -  iV  ('l"i.,'  =  (r-'(i)S  -ik-  ^^^u]\ 

o  /  "^  o  • 


arU,  Y~   U 

IZ  _  J^ £  arP 

1_U  1-U        '^^'^Ir 


rANO 


,AsO 


-1 


-lU,       =    (#)°    U,    +    irJn    C    C,    -    (C^    +   U    (1-U    ))D    -^arU, 
1    In         ^D   U      1         ^D'C      o   1  o  o  o        0  Ir 


-1    P   -1 

-    Y  C^D   -^    arP,     , 
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The  form  of  these  equations  suggests  changing  to  new 
dependent  variables 


Vi(./,r)  =  - 


-1 


(27)       V2(n,r)  =  y"  Pi(*(,r)  +  V-^  (,|,r  )U^  ( ,0 

Equations  (26)  then  reduce  to 

p  . 


aU 


'iV 


1»-|   1-IJ 


XV, 


/     2G'-  \         C^ 

-  1  +  w-^\   rV,   +  -^^  rV„   + 
U  D  I   Ir   „2^  2r 


o  o 


U  D 
o  o 


2U  -1  A     1-U 


o   '"o  ,    'o  /AsO    Y"l   o  /A^o 


(28)   nV 


U        O 

•-  o 

aU   ( 


+ 


(-)' 


aU   f 


U    'D'U    2   U   ^D'C 
o  o 


-1 


^D^C  y2   2   '-^31 
o        -' 


l^r/  i 


rV.  ,  -  rV^^  +  XV.  [3a-2Y"  (1-a)] 


2r 


'iV 


( 


o   -  ^-Tt-  S  ^^1       -  rV^   +  XV, 


n1 


3a  -  — Y  (1-a) 
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We  introduce  the  new  independent  variables 


(29) 


[ ^Tnrr-T  ^  -  log  r  =  y  -  1 


-Og  r 


In  which  H„   is  the  value  of  ri     for  which   1-U   +  C 


0, 


o     o 
that  is,  D  (Ii„  )  =  0.   iVote  that  the  lines   x  =  constant   are 

the  particle  paths  of  Guderley's  flov;,  see  (17),  and  that 
y  -  0   is  a  Mach  line,  see  (19).   Increasing   t,  which  corres- 
ponded to  increasing  y^  ,    now  corresponds  to  decreasing  y. 
With  these  variables  equations  (28)  become 
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2         ? 
2C        C 

o  ly    U    Ix    tt2   2x    1  o     30 

o 


(30)       V^y  +  V^^  =  -tx^V, 


3y    Ix    ^3  1 


where 


i2TJ  -1      1-U      „  .      ,  C 


U  o  o  ) 

{  o 

(31) 

O 

^3       Y-1     Y~l 

These  are  the  differential  equations  that  are  to  be 
solved. 

5«   Boundary  Conditions  at  the  Shock, 

We  assume  the  path  of  the  shock  is  represented  by 

M  =  H^(l  +  e  H^  (r)  +  ...) 

(32)       ^    r 

or   t  =  r  K^(l  +  eH^(r)  +  . . . )   , 

The  piston  starts  to  move  at  time   t   +  et,   from  the 

o     1 

position   r  =  R   +  eR-,  .   VJe  assume  that  the  shock  issues  from 
the  piston  at  the  same  time,  that  is, 

(33)    "o  +  et^   =    (R^  +  eR^)^  K^d  +  eH^  (R^  +  e\)    +  ...). 
On  taking  first  order  terms  in   e;  ,  v;e  obtain 

H,  (R  )  =  — aR.  R  ^~^ 

1   °    H  R      1   ° 
o  o 


1^. 


or  using  the  zero  order  result,  t   =  R   H  , 
°  •*   o     o    o' 


(31+) 


H.  (R  )  =  R 
1   o      o 


X-1 


t     aR, 

t     R 
o     o 


In  this  relation,  see  Section  l\.,    t,   and   R,   can  be  chosen 
arbitrarily. 

Along  the  shock  path  vie  have 


(35) 


dr 
dt 


r  / 

=   z   =   a  TT  Z      or   Z  =  Ad  log  r/d  log  t 


and  hence,  by  (32), 


(36) 


1 


acrH 


1  + 


Ir 


earH^   + 
Ir 


1  +  eH 


1 


where   H-j^(r)  =  dtl-^/dr.   Substituting  (36)  in  (12)  and  us  in^ 
(15)  and  (20)  we  find  on  the  shock, 


(37) 


U  =  U(H^(l  +  eK-^(r),r)  =  (1-iJ.'^)  (1-earH^^.  . .  ) 
C  =  C(H^(l+eH^(r),r)  =  M-J^I+^l^  (1-scrH^^^. .  .  ) 
P  -  F(H^(l+eK-^(r),r)  =  (l-iJ,^  )p^  (l-2earE^^.  . .  ) 


Cl.  =  ^   ,  see  (7 ) , 
C 
Expanding   U,  C,  P,  in  powers  of   e   as  functions  of   e 

for  fixed    n,  r   we  find  using  {Ik) 

U(H^(1+En^(r)),r,£)  =  U^(H^)  +  eU^(H^,r)  +  eH^(r)[|^j 


ii 


(38) 


UJHJ  +  eU^(K^,r)  +  ei^ 

o 


o'-O'     -"l^'-o'^  '   •  "'yD 


H-^(r)... 


-i=H 


C(H^(l+eH^(r)),r,e)=C^(H^) 


B 


l+eC^(Ii^,r)+e  ^    H^(r)... 

I   o 


P(H^(l+eH^(r)),r,£)  =  P^CH^) 


15. 

(E  \ 

/   o 


By  comparing  (37)  and  (3j),  and  using  (15),  we  obtain  the 
boundary  conditions  for  U-,  ,  C-,  ,  P-,   at  the  undisturbed  shock 


curve   n  =  H 


1 


Ui(H^,r) 


(39) 


^  o' 
Oi(H,.r)  =  .^^ 

^       O  ' 


=K 


H^(r)  -  (l-n.'')arH^^(r) 


H^(r)  -  (l-ii'^)arH-L^(r) 


»;=H. 
(   0 


p  fH   r )  =  -*  — 
^l^^o'^^     \D  1 
\  o  ■ 


H-^(r)  -  2arH^^(r). 


f;=-tl 


IVe  suppose  now  that  the  undisturbed  shock  curve   •/  =  H 

corresponds  to  j   =   y(K  )  =  Y  ,   The  three  shocli  conditions 

'^    ■•   o     o 

(39)  expressed  in  teims  of  the  variables   V-,,  \',    V„ ,  x   and 
y,  see  (27)  and  (29),  are  then 


(1+0) 


V^(x,Y^)  =  p^II(x)  +  q3_H^(x) 
V2(x,Y^)  =  P2H(x)  +  q2\(^) 
V^(x,Y^)  =  p^H(x) 


vjhere 


(1+1) 


ik2) 


H(x)  =  U-^ir) 

,.  2^-l  -2  /  ^c 


o  - . 


/   o 


^o^ 


/  -1  ^o  ,    . 
P2  -   -"^   —  ^  —  D 

I  O     [.L     O 


/M=H_ 


Qr:.!    /; 


*  ;  .' 


,  i/.    '    y 


«  -.-^1. 


16. 


„   -  f   _2 0  ,     1 o\ 

I   ( 


q,  =   -a/[i' 


'1 


q^  =  (2y"   -  M-"   +  l)a. 


These  boundary  conditions  contain  an  umdetermined  function 
H(x)   wliich  can  h.o^^^ever  be  eliminated  in  an  obvious  way.   The 
equations  which  result  can  be  written  In  the  form: 


{k3) 


p^q2V^(x,Y^)  -  p^q^V2(x,Y^)  +  (p^q^ 


Pp^a^l^^'^'o^  -  PiP3V2(x,Y^)  +  (P^q2  - 


with  the  additional  condition  V-,  (X  ,Y  , 

3   o'  o 


by  {3k),    (^0)  and  (kS). 


p^q2)V^(x,Y^)  =  0 


P2qi)V3^(x,Y^)  =   0 

^  A-1   "^1      ^1 
P3^o      r  ■  ^  ~ 
L^o 


6 •   Boundary  Conditions  at  the  Fiston» 

We  assume  that  the  piston  velocity  is  given  through 
equations  (21)  by 


ihh) 


u  =  u  (,, 
p  1 


U^(n)  +  eU,^(.i). 


■IP^'/ 


The  connection  between   U  ,  U,  ,  R    and  R,   will 

o '   Ip '   o        1 

flven  later.   However  from  (23)  we  have 


be 


U(.,j,R^(>P  -H  eR3_(,^),  e>    =  U^(.j)  ^  eR^(:,)  ^(U^(.,))  ^  sTJ^(.pR^: 

=  U^(.)  -.  sU^(,,,R^(.|)) 


or  using  {l\M-) 


.     I 


fill,.' 


17. 


(1+6)      U^(.,,R^O-|))  =  %(')). 

In  the  (x,y)-plane  the  undisturbed  piston  path   r  =  -^n^*^/^   ■^^ 
given  according  to  (22)  and  (29)  by  the  line   x  =  constant,  say. 


ik7 


X  =  X  . 
o 


Expressing  (li-6)  in  terras  of  the  variables   x,  y,  and  V,  by 
(27)  and  (29)  we  have 


(48) 


U.     (>](y)) 
^l^^^o'^^^^    =  -   Ull-U    )      • 


We 


next    determine      U,     ('-i)      in   terms    of      U    (n),    R    (n) 


-X. 


and  R,  (;/,r).   From   r  =  R  (r  't ) ,  see  (2),  we  have 


dr 

d 


dR 


t    r  ds   ^  ^  dt    t^  ' 


V,         -uu-  ■     •     TT   ^t  dr 
hence  we  obtain  using   U  =  —  -jr   ; 


(49) 


r  dt 
dR  aR  ('.,)  U  ( 
d'  "  '" 


1-U  ((,}  • 
P  / 


From  r  =  Rp(.|)  =   R^(/^)  +  eR^(r,),  U  =  Up(i^)  =  U^(>;)  +  eU^p(K^) 
we  obtain  to  first  order 


dR 
c 


+  e 


dR^  ^  a(R^(»))+eR^(t;))  I    ^ 


-^'^V^ 


aR   U 

o   o      ae 

o    I  ^    o 


7 


R  U 


U^(r,)+sU^p(.^)_ 


U^R,  +  Y^Y^  U,   +  U  R  TT 
o  1    1~U   Ip     o  o  Ip 


or 


dR.    aU 


U. 


ITT  ^1  ^ 


iE. 


(i-u^) 


■o  aR 
2    o 


1 


■ir.-rr.-T-:     !"■ 


i:  I  -'  -  1  '^  .  X 


ML  » 


;■''!'. 


18. 


or 


(50) 


u. 


ip^l' 


R 


2r 


dR  dR 

I)  x-^  +  RVi    :r^ 
I  d\j  1(     drj 


7»   Conditions  at  the  Junction  of  Shock  and  Piston, 

From  (31;)  and  ([{.l)  vje  have  the  value  of  H(x)   for  x 


o' 


(51) 


H(X  )  =  R  ^"^ 
o     o 


t       R 
o      o 


From  the  first  shock  condition  (1|0)  and  equations  (21),  (1|8 ) 

and  (^0)  vie   obtain  a  relation  between   K  (X  )   and   H(X  ), 

X  ^  o  ^  o  ' 


P3_H(X^)  ^  q^H^(Xj 


Mi-iyH^ 


-1  ^  /dR 

H  [   1      +  R  H  (-^ 
o\  d>|  /_,-     1  o  \d'i 


/dR^ 


/  'H 


or  using  (22)  and  (29) 

(52)       P.IUX  )  +  q  H  (X  )  =  (-r^l       -  R, 
1    o'    ^1  x'  o     \d7  J       =Y  ^ 


'dR- 


where  we  noI^r  consider   R-,   as  a  function  of  y   instead  of  n  , 
see  (29). 


o-" 


8,   The  function   G(x,y). 

The  problem  has  now  been  reduced  to  the  following:   the 

three  functions   V-,  ,  Vp   and   V    satisfy  three  homogeneous 

1   ^        3 

linear  partial  differential  equations  (30)  in  x   and   y 
whose  coefficients  depend  on  only  one  variable   y   and  are 
analytic.   One  boundary  condition  (i|8)  is  given  on  the  piston 
path  X  =  X   and  two  conditions  (Ij-3 )  on  the  shock  path  y  =  Y  , 
Both  these  lines  are  tiraellke,   V/e  may  expect  that  there  exists 
a  unique  analytic  solution  of  this  problem  if  the  boundary 


( j;s , 


'■  r 


19. 

conditions  are  analytic.   Once  its  solution  is  known  the 
function  H(x)   can  be  foiond  from  the  last  equation  of  (i]-0). 

The  form  of  equations  (30)  leads  us  to  introduce  the 
function   G(x,y)   such  that, 

(p3)  V   =  G 

1   y 

iSk)  G(x,Y^)  =   0. 

Then  we  may  integrate  the  last  two  equations  of  (30)  and 
obtain,  by  (5/4.), 

(55)  V2(x,y)  =  V2(x,Y^)  -  G^(x,y)  +  [i^GU,j) 

(56)  y^U,j)    =   V^(x,Y^)  -  G^(x,y)  +  ti^G(x,y). 

Substituting  (53 )>  i5S)    and  (56)  into  the  first  equation 
of  (30)  and  using  the  shock  conditions  (I(-0)  we  find 
2C         C  /         C 

o  \  o  / 

-  P3LoH(^)  -  P2  ^  \(^)  -  ^-2  ^  ^xx^^^- 

o  o 

In  this  equation,  v;e  recall  D  ,  U  ,  C  ,  K   and   L    are  all 
'■  '  O'   O'    o'   o  o 

functions  of  y.   On  multiplying  through  by  the  quantity  ^ 

o 

v/hich  has    a   finite    value    at      y   =   0,    D      =   0,    we    obtain   finally 
(58)      yG        +  a(y)G      +  b(y)G        +  c(y)G   +  d(y)G      +  e(y)G 

=     i:(y)II(x)    +  m(y)H    (x)    +n(y)H      (x). 
With 


!    ^   '  1 


>"  \  1 1 


20. 


a(y)  =  g-  K 

b(y)  =-^,r^ 

o   o 

c(y)  =  §-  ti3  L^ 

d(y)  =  -  ^  L^  -  tx^  -§ 


(59)  "o  \  "    "  U 

C2 


o  U 

0 

i(y)  =  -  g-  P3  Lo 

C2 

o    U 
o 

n(y)  =  -  ^  ^2  tI 
o    U 
o 


None  of  the  functions   a(y),  b(y),  c(y)  ...   etc,  is  singular 

for  V  <  Y  . 

-  -  o 

The  boimdary  conditions  on   G   for   y  =  Y   are,  by  (53) 

and  (5i|)  using  (1^0), 

(60)  G(x,Y^)  =  0,   C-^^(x,Y^)  =  p-j_n(x)  +  qiH^(x) 
and  for  x  =  X   from  (^3),  (5'l-)  and  (i|8) 

(61)  G(X^,y)  =  -  '   TT  n-?T  \   dy. 


^o   °    °' 


U  (1-U  j 
o     o 

Two   conditions   are   iraposed   on     H(x)      through    (5l)    and 


(52), 


■i•i^  iia,r,; 


21. 


:62) 


H(X  )  =  R 
^  o     o 


X-1 


ll  _  ^\ 


t 


R 


p_H  (X  )  +  q^H  (X  ) 

^1X0      ^1X0 


.-dR^ 


y=Y^ 


-  ^1- 


Nox-7  the  problem  is  reduced  to  solving  a  second  order 
hyperbolic  differential  equation  (58)  for   G  as  a  function 
of  X   and   y.   It  involves  in  the  non-homogeneous  term  an 
unknown  function  H(x),   On  the  line  y  =  Y  ,  the  two  con 
ditions  (60)  for   G   and   G   are  imposed  and  one  condition, 
(61)  for   G,  on  the  line   x  =  X  ,   The  quantity  H   and  its 
first  derivative  are  known  at  the  point   x  =  X  ,  equations 
(51)  and  {S2). 

An  indication  that  these  data  determine   G(x,y)   and 
K(x)   is  given  by  the  follovjing:   If  the  "piston  function" 
G(X  ,y)   can  be  represented  by  a  power  series  around  y  =  Y  , 
then  it  i;,  possible  to  determine  a  power  series  expansion  for 
G(x,y)   about   X^,  Y^. 

V7e  wish  to  determine  the  behavior  of   G  at  the  c  enter, 

X  =  CD  ,   We  would  like,  in  fact,  to  determine   G(x,y)   in  the 

domain  X  <  x  <  co ,   Y   <  y  <  Y   where   y  =  Y    corresponds 
o-       •'oo-'^-o         "^     oo        -^ 

to  the  position  of  Guderley's  reflected  shock  at  n  =  K   • 
Beyond  y  =  Y    v;e  would  need  to  fit  a  flow  satisfying  appro- 
priate shock  conditions  at   y  =  Y   ,   This  aspect  of  the  prob- 
lem is  not  discussed  here  since  we  are  only  concerned  with 
finding  out  if  the  perturbation  approximation  breaks  down  as 
the  flow  approaches  the  center.. 


22. 

We  examine  now  the  flow  picture  in  the  (x,y)-plane.   The 
differential  equations  of  the  characteristics  are  found  most 
easily  from  equation  (l-(-2): 


^dy;    U^   dy/-  .2 


0 


or 


C  C 

df  "  D~§"  ^'^o  t   (i-Uq)]  =  u  Lc  4i-u  )J  • 

•^  o  o  0   o-     o 


Thus  the  line   y  =  0  vjhlch  corresponds  to   (1-U   +C)=0 
•^  ^  ^      o    o 

is  a  characteristic. 


:.-;■    ^k:.;- 


?n.:.i.! 


J-        /   .-  - 


C. 


n 


23. 


Tlrae 


Piston 


(^o'^o) 


(Xo,y;L) 


(x^,o) 


(Xo^Y^o^ 


Characteristics 


Particle  paths 


Figure  1. 


'V  !■:*;;.! 


i 


\ 


\l 


-'i.,.Y,,_,. 


•1  .J  .•*.1  J  i  ■.  1' 


2l+. 

From  Figure  1  we  see  that  since   y  =  0   is  a  character- 
istic, the  value  of   G  at  any  point   P,  with  X  <  x  <  oo, 
0  <  y  <  Y  ,  depends  on  the  boundary  values  only  in  the  range 

X^  <  X  <  X-,   on  Y  =  Y   and   y,  <y<Y   on  x=X   where 
0--1      -     o       •'l-''-o  o 

y-,  >  0,   Hence   H(x)   for  all  values  of  x   is  determined  by 
the  piston  motion  up  to   y  =  0,   From  the  figure  it  is  clear 
that  any  slight  discontinuity  in  the  piston  motion  near  y  =  0 
might  produce  large  effects  for  large  values  of  x.   '-'e  there- 
fore assume  that   G(X  ^y)   is  analytic  for  y  <  Y   and  expect 
that  a  unique  analytic  solution  exists.   For  the  time  being 
we  consider  only  the  segment   0  <  y  <  Y  ,   For   y  <  0  we 
simply  continue  the  flow  analytically. 

Part  II,   Application  of  the  Laplace  Transform 
9 •   The  Transformed  Problem , 

In  the  partial  differential  equation  (58)  for   G(x,y) 
the  coefficients  depend  on  y   only  and  not  on  x.   This 
suggests  introducing  the  Laplace  transform.   The  use  of  the 
Laplace  transform  is  usually  restricted  to  cases  in  which 

the  normal  derivative  at   x  =  X   is  eiven  as  well  as  the 

o 

function  itself.   However^  the  property  that  the  lines   x  =  X 

and  y  =  Y   are  both  timelike  and  that  the  equations  are 

hyperbolic  can  be  used  to  derive  a  relation  between  the 

normal  derivative   G  (X  ,y)   and  the  given  boundary  condi-r- 

yi      o 

tions . 

¥e  assume  a  unique  analytic  solution   G(x,y),  H(x)   of 
equations  (58-62),  Part  I,  exists  such  that 


-n 


2^, 

(1)  g(y,k)  =  J   e-^^^'-^o^  G(x,y)dx 

X 
o 

r°°  -^(x-x  ) 

(2)  h(k)  =  J    e       °   H(x)  dx 

exist  and  converge  absolutely  for  Re  k  >  K'        say. 

If   G(x,y),  H(x)   satisfy  (58-62),  Part  I,  then  the 
transforms   gCyjk),  h(k)   satisfy 

"  I  p 

(3)  yg   +  (a+bk)g   +  (c+dk+ek  }g 

=  (i  +mk+nk^)h(k)  +  P(y)  +  kQ(y) 


X'lith 


a(y)  =^K^ 
o 


2 


b(y)  =   -   ^jr^ 


D   U 
o   o 


°  2 

d(y)  =-i-L-   1^2  Ji 
o  V        iJ 


C2 

o  IJ 
o 


i(7)    =  -  f  P3  Lo 

0   -^ 

m(y)  ^  -  T?-  Pp  -§ 

^o  "^   U"^ 

o 

n(y)  ==  -  ^  q2  -f  . 


o 


U' 


-h' 


--■  -1 


a/i'V  '.: 


I  '_-, 


■\     I 


:'-    r:,i-    1; 


26. 
Here 


ik)      P(y)  =  bGy(X^,y)  +  dG(X^,y)  +  eG^(X^,y)  -  mH(X^)  -  nK^(X^) 

(5)   Q(y)  =  eG(X^,y)  -  nH(X^), 

G  (X  ,y),  G(X  ,y)   are  given  by  (60 )  and  (61 ),  Part  I,  and 
H(X  )   and  H  (X  )   by  (62)  while   G  (X  ,y)   will  eventually 
be  determined.   We  note  here 


(6)  ^^(^o^^'o^  =  ° 

by  (60)^ and  (6l),  Part  I, 

Prom  (1)  and  (2)  and  the  boundary  conditions  (61 )  vie 
obtain 

(7)  g(Y^,k)  =  0 

(8)  gy(Y^^k)  =    (p^  +  kq-^)h(k)  -  q^H(X^). 

Since   G(x5y)   is  regular   in  y   we  also  have   g(y,k) 
analytic  in  y.   Equation  (3)  has  exactly  one  singular  point 
y  =  O5  and  thus   g(y,k)   is  regular   for  all   y  <  Y^   if 

(9)  g  is  regular  at  y  =  0. 

Since   g   is  a  solution  of  (3)  it  can  be  vjritten  with 
appropriate  functions   A(k)   and  B(k)   as 

(10)  g(y,k)  =  A(k)g-j_(y,k)  +  B(k)g2(y,k) 

+  h(k)g^(y,k)  +  g|^(y,k) 
in  which   g-,  (y^k),  gpCy^k)   are  two  independent  solutions  of 


\     ^. 


:    ■.      I  7.     A^     r 


ix,.^)^^ 


j':-n    J* 


.-.,■    f  r\ 


:iT  ■ 


/     ■'    V  ■- 


'ICj. 'i-j::*:; 


I  ' 


-1   ,  ' 
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the  homogeneous  equation 

(11)  X(g)  =^  yg   +  (a+bk)g'  +  (c+dk+ek^)g  =  0 
while   gT(y,k)   is  any  regular  solution  of  the  equation 

(12)  /-.(g)  =  i  +  mk  +  nk^ 

and  gi  (y^k)   is  any  regular  solution  of 

(13)  ^'(g)  =  P  +  kQ. 

We  specify  no\-j   which  solutions  vxe  consider  for   g-,  ,  gp, 
g„   and   g,  ,   The  indlcial  equation  of  (11)  for  y  =  0   is 

5(5-1)+  (a   +bk)5=0 

where   a   =  a(o),  b   =  b(o).   Equation  (11)  therefore  has  one 

solution  that  is  regular  at   y  =  0,  and  one  solution  v/hich 

-a  -b  k+1 
behaves  like   y  as   y  — >    0   except  when   -a  -b  k+1 

is  an  integer. 

Here   g, (y,k)   is  the  regular  solution  of  (11)  with 

(1^)  &-^{0,]z)    =   1 

and   g  (y,]')   is  the  singular  solution  with 

-a  -b  k+1 

(15)  SpCy.k)  =  y    °    °      (1  +  o(y)) 

near   y  =  0,   Then  g^  (y^k)   is  defined  except  when   a  +b  k 

is  a  negative  integer  while   gp(y,k)   is  defined  except  when 

a  +b  k   is  a  positive  integer, 
o   o        ^  " 

¥e  turn  next  to  the  solutions  of  (12),   Since  the  non- 


28. 


homogeneous  term  Is  regular  In  j      the  regular  solutions  of 
(12)  form  a  one-parameter  family,   VJe  fix   g  (y,k)   for  dif- 
ferent ranges  of   k  by  the  additional  conditions 


(16) 


^(0,k)  =  T-,^   for   Re  k  >  0 


;  (Y^,k)  =  -Up    for   Re  k  <  0, 


From  general  theory  the  function  g^{j,k)      is  defined  except 

Xrfhen   a   +  b  k   is  a  negative  inte-'^er, 
o     o  '^^ 

Similarly  we  might  fix   gi  (y,k)   by  the  condition 

g^(0,k)  =  T^i^     for   Re  k  >  0 

(17) 

gj^_(YQ,k)  =  i^j^^_  for   Re  k  <  0. 

Thus   g-[_(ysk),  g_(y,k),  gL(y,k)   are  all  regular  in  y 
but   gp(y5k)   singular  at   y  =  0,  provided   a   +  b  k   is  not 
an  integer.   Since   g   is  analytic  in  y^  we  have 

(18)  B(k)  =.  0. 

From    (7)    and    (8)    \<ie    can   determine    the    other   coefficients    in 
(10). 

A    {k) 

(19)  A(k)    =  ^^^7j- 


(20)  h(k)    -      ^ 


h,  (k) 


D(kT 


v/here 


;21)      A-j_(k)    =   g^g^y   -    g3g^y   -    ci^i}{Xjg^    -    (p^    +   kq^)r^ 


29. 

(22)  h^(k)  ^  i^i^^^  -  i^i^^  +  qiH(X^)il 

(23)  D(k)  =  i3iiy  -  iig3y  +  (Pi  +  kq^)i-j_  . 

Here  the  bar  means  evaluated  at   y  =  Y  .   Note  that   D(k)   Is 
Independent  of  the  value  of   go  at   y  =  0,  see  (16). 
Thus  V7e  have 

A^(k)  h  (k) 

(2/i)  g(y,k)  =  p(,^)   gl(y,k)  +  5^]^  g3(y,k)  +  g^(y,k). 


Since  g^ij 3^'-)      drops  out   g(y5k)   Is  defined  for  all  value 


s 


of   k  unless   a   +  b  k   Is  a  negative  integer  or  if   D(k)  =  0, 

00  °  \     I  • 

At  this  stage,  the  transform  g(y5k)   cannot  be  deter- 
mined directly  from  the  boundary  conditions  given  for 
G(Xjy),   For,  the  function   g,  (y,k)   depends  on   P(y),  see 
(13)5  which  in  turn  depends  through  ().j.)  on  the  function 
G  (X  .y).   However,  from  the  behavior  of   g  for   Re  k   large 
and  positive  vie  can  determine  a  relationship  between  G(X  ,y) 
and   G  (X  ,y).   It  will  then  be  possible  to  determine  some 
properties  of   G(x,y)   from,  the  properties  of   g(y,k). 


9 .   As^/mptotic  Behavior  of   g(y,k) , 
We  have ,  see  (1 ) , 


pCO 

,10  -\ 


-k(x-X^) 
e 

X 


g(y,k)  =  f    e      °   G(x,y)  dx. 


o 

For   Re  k   large  and  positive  the  main  contribution  to  the 
integrand  of  this  expression  comes  frora  x  =  X  .   Making  use 
of  this  fact,  ^^^e  find  for   Re  k  — ^  +  00, 


>  ,-,f^- 


30. 

(25)  g(y,k)-  i  (^(^0'^^  ^   °(^0  • 

If  we  use  expression  {2l\.)    for   g(y,k)   we  should  find 
the  same  asymptotic  behavior  for   g(y,k).   With  this  goal  in 
mind  we  determine  the  asymptotic  behavior  of   g-,(y,k), 
g_(y,k),  g|  (y^k).   This  is  done  purely  formally  here  and  the 
results  are  proved  in  [L[.]  , 

Using  the  usual  procedure  for  finding  asyinptotic  ex- 
pansions we  set 

(26)  g^  =  e"^  ^J'l(y,k) 

and  substitute  in  the  differential  equation   ^-Kg-,  )  =  0, 
see  (11),   VJe  obtain 

i;  ,1  It      , 

(27)  Ji<-^   +  (a  +  bk  -  2w  yk)\l/-^^  +  (c  +  (d-yw   -  aw  )k 

+  (e  -  bw   +  y(w'  )^)k^)\|/^  =  0, 

The  function  w(y)   is  determ.lned  by  setting  the  coefficient 

2 
of   k    equal  to  zero, 

(28)  yw'2  -  b(y)w'  +  e(y)  -  0, 

Since   g-,(y,k)   is  regular  at   y  =  0,  w(y)   is  chosen 
regular,  so  that  from  (28)  we  have 


(PQ]  dw  _  b(y)  -   Kb"(y)  -  h.o(y): 

dy  2y 


Then  the  differential  equation  for   \1/-,  (y,k)   is 


30)      y\|/^  +  (a  +  b^k)\l/^  +  (c  +  d-j_k)\l;^  = 


0 


'to 


.-.";'■•.,■.•      ■"!      .' 


/ 


V;  .  I. 


31. 
where 


2C 


b-,  (y)  =  Kb^(y)  -  k-e{j)j   =  ^  T^  (1  -  U  ) 


31) 


o   o 

t         " 

d-^Cy)  =  d(y)  -  a(y)w  (y)  -  yw  (y)i 

From   g-j^COjk)  =  1,  see  (li;),  and  (26)  we  have 
w(0)  =  0 


>^^(0,k)  =  1 


°   °  by  Part  I  (59) 


and  therefore 


-y 


(32)  w(y)  =  j   bCr)  -H'^in-keU),   ^,/. 

"0 

Since   b(y)  >  0,  b  (y)  -  I|-e(y)y  >  0,  e(y)  >  0,  (see 

Part  I  (59)),  the  function  w  (y)   is  always  positive  and 

different  from  zero  in  the  range   0  <  y  <  Y  ,   At   y  =  0, 
,       e 
-T-r  =  -  r—  <  0.   The  function  w(y)   is  therefore  monotonia 

o  2 

increasing.   The  condition  b  (y)  -  i|e(y)y  >  0   is  exactly 

the  condition  that  equations  (53)  for   G(x,y)   is  hyperbolic. 

The  condition  b(y)  >  0   is  connected  with  the  fact  that  the 

two  lines   x  =  X   and  y  =  Y   are  both  timelike.   As  we 

o       "^     o 

shall  see  later,  the  possibility  of  solving  the  problem  by 
using  the  Laplace  transform,  rests  on  the  fact  that  w(y)   is 

m.onotonlc  Increasing, 

2 

From  equation  (30),  v/hich  contains  no  term  in  k  ,  an 

asymptotic  expansion  for   \|/-,   is  found  by  setting 

^nT(y)    ^np(y) 

(33)  ^i(y,k)  -^^-Lo(y)  + -nz —  + ^~-  +  ..,     . 


(  -s 


r  .   a"  •  J  I" 


32. 


Differential  equations  for  ^-,q{j)}    ^-|-|(y)j  etc,  are 
found  by  direct  substitution: 


^1^10  +  ^1^10  =   ° 


'3k) 


The 


t  'J         I 

l^^lii    l^ln    -^^l^n-l     ^l,n-l     ^l,n-l 

condition  (31)  yields  initial  conditions  for   \}''   (y) 


(35)  ^10^°^  =  ^ 

t^^(O)  =0     n  >  1. 

Equations  (3^)  and  (35)  are  easily  solved. 
A  formal  expansion  for  g-yij,^),    see    (12),  can  be  found 
in  a  similar  way.   We  set 

il'.T(y)   ■■■ 

(36)  g^(y,k)  ^  ^^Qij)    +  -^3^ +  ...  -  A^(k)g^(y,k) 

where   \l/„„(y),  .,,,  H'o^(y)   are  determined  by  straight  sub- 
sbitution  in  (12),   '■'-■'(go)  ~  -i-   "^  ™k  +  nk  ,    and  A-,(k)   is 
left  undetermined,   V."e  find 


'37)  +   e(y)\l/^Q  =  n 


m 


+b(y)4'^Q  +  d(y)i^^Q  +  e{j)^^-^ 

y^30  +  a(y)vlf^Q  +  b(y)\l/^-^  +  c(y)i^^Q  +  d(y)^^-^  +  e(y)\!/^2  "  ^ 

^^(y)^3,n-i  ■"  ^(y^^,n  =  °- 

Since   e(y)  =  -  "§ p   ^^®®  Part  I,  (59))  does  not 

°  ^o 
vanish  on   0  <  y  <  Y^   the  functions  "^-iQ^y)  >    ^-.j^ij),    etc. 


1 


-C 


;.-    n 


ri2 


{-;...'/ 


'  -'>.-.  'I' 


I.    ;     .,.:- J•,':.^.:;• 


iT 


..    «    «  «         j 


1  1 


O,'^''    '' 


■I- 

i   -  ■  ;  . . 
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are  all  determined.   We  now  fix  A    (or  t    and  t„   in 
(16))  asyrflpt otic ally  by  requiring 

g  (0,k)  -v.  0     for  Re  k  >  0 

(33) 

^|/..  (Y  )    ■ 
g^(Y^,k)  r^i^^^iY^)    +     ^^^  ""■     +  ...)   for  Re  k  <  0. 

In  a  similar  way  we  fix   g,  (y,k) 

(39)      gl4.(y'^)  '-  k  \\o^y^  ""    k ■"•••/     ^""^  Re  k  >  0 

where 

(^0)  e(y)i['^Q  =  Q(y) 

()+i.i)       b(y)\o  "^  ^('^^\o  "^  e(y)\i  =  P(y) 

(J+1.2)    j^'^Q   +   a(y)\l/^Q  +  b(y)\l;^^  +  c(y)\|/^Q  +  d(y}\l/^^ 

+  e(y)\l/^2  =  ° 

Ca.n)     y^i!,n.2  ■*■  ^^y^\,n-2  "^  ^(^^)\,n-l  "^  =^y^\,n-2 

-^  'i(^^^\,n-l  -^  ^(y^\,n  =  °- 

For  Re  k  <  0,  g,  (y,k)   is  fixed  by  a  formula  analogous 
to  (38).   Combining  (26),  (33),  (36),  (39)  we  find 

T.(y,k) 
(i|2)     g3(y,k)  -  4'3o(y)  +  -^^^ 

T,,  (y,k) 


1  /      ^1,  ^y,i<^  >* 


(■  c  ■' 


Ti       r-. 


c  '■ 


3i  Si'i:    ••ui'i 


-to. 


7      ■.^■ 


r-A       :al 


^  •  \ 


>i 


(y:..,'  - 


•;'  \   r  , 


)  r    •>.■    r.  4^:'Z)  0    -^     f   jv  (  "i")''''   -^    r 


•'  ■  ^    :!• 


i  ^..<  "i  r. 


I  -■■' 


'oi-jnoijsn^    .P  Uf!.;"!v'l    ^    ■•;;-!    i?3;^:_    g. 


3i;. 


where 


•  •  • 


([|3)     T^(y,k)  =  *3i(y)  +  -^^^ —  +  ... 

T^(y,k)  =  ^^^(y)  +  ^ +  ... 

§re  all   0(1)   with  respect  to  k.   It  has  been  shown  in  [[|.] 

that   T^,  T   and  T,   are  all  bounded  for   |k|   sufficiently 

larce  and  Re  k  >  k   where   k   can  be  chosen  arbitrarily 
°  p  p  "^ 

but  cannot  be   -co  , 

Substituting  (Ij.2)  in  (19-21)  we  find  the  asymptotic 

expansions  of  the  coefficients   A-,(k),  h,  (k),  D(k) 

D(k)  =  i3iiy  -  i3yi;^  +  (Pl+kqi)i-L-^  k(w'(Y^)^3Q(Y^) 

-kw(Y^)  -kw(y^) 

+  qi)e         ^lo^^o^  ^  ^        T^^Y^'k) 

{hh)    h^(k)  =  i3_i^y  -  i^i^y  +  qiH(X,)ii  -  (w'(Yj^^o(Y^) 

-kw(Y^)  -.      -kw(Y^) 

-^  q^H(Xj)e         ^^^(Y^)  +^  e        T^(Y^,k) 

Ai(k)    =  i^i3y   -  i3i^^y  -   qiH(Xjg3    -    (P3_+kq^)i^ 

^-q-LH(Xj^30   -    q^^^o    -^iT^(Y^,k) 


where  T^^,  T,,  T^   are  all  bounded  in  the  same  sense  as   T^^, 
nd   T^. 
We  note  here  in  particular  the  expression  for  T„(Y^,k) 


T.   and   T,  . 
3        W 


("•^  ,. 'i^ 


•  •  •       '' 


.'fiiod   'tiisrl   vtl 


T 


.  OT!- 


[y.)C    .^A)^ 


'io 


(o" 


I      '.:  i  .'1 


<     -^      ( 


9     -r      i 


vc 


y}  w) 


1 - 


o  '  ' 


•'fh  ■ 


s:    ( 
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ikS)   T^(Yo,k) 


^1,1  (Y^) 


=  l(y<^o)  ^^TT^*-)  (♦30(^0' 


+ 

00 


\l(^o) 


n-1 


^ZI;;^   H  ^V^3,n-r  -  V^3,n.r 


n=0 


0 


-  Pl\,n.l  -  ^l\,n  -  ^l^(^o^*3,n^ 


where  the  bar  indicates  evaluation  at   y  =  Y  • 

From  {l\l\.)    and  {2k)   we  obtain  the  asymptotic  behavior  of 
g(y,k)   for 


T^{Y^,k) 


(I1.6) 


g(y,k)'N.  -;3-^^ V(Y^,k) —  k 

-     -._  ^e^^o'^^ 


^10  -^ 


T-L(y,k)" 


k(w-w) 


^10^-^  ^30''^l^ 


T^(Y  ,k}   k   ^30'*' 


T,(y,k)" 


^l 


\o  "■ 


\(y,k)' 


K 


T^(Y^,k) 


^10^-qi^(-^o^^3Q"^l\o^      k     ^15  k(w-w) 


^10^""  ^30  ^  ^^  ^ 


1+0  ( 


^1' 


^'] 


,  ho(^'\o  -^  qi^^^o^  .^. 


^10^-^  ^30  ""  ^1^ 


30 


H(^^]  ^^\o(i^°(|^ 


I     "^ 


-  . » .  + 


.<.-»•   — 


M 


.• —  ^ 


»'^  r»  '"Y        >     *    '  i    I  . 


■T~n, 


CO 


/^i^  - 


^■x^ 


'io   'xoiViirlec    S^.' 


V  '• 


n 


X"     u;. 


-  i'  I 


L     ^ 


■:;.;!- 


Ji 


36. 

From  (l|-0),  (36)  and  (5)  we  have 
ikl)  \oij)    =   Q(y)/e(y)  =  G(X^,y)  -  ^   H(X^) 


=  G(X^,y)  -  \l/^Q(y)  H(X^). 


Using  also  G(X  ,Y   )    =   0,    (6),  the  asymptotic  expan- 


sion for  g(y,k)   reduces  to 

T  (Y  ,k)^i/    k(w-w)  r      -, 
( -w    Wi^  _+n-  )     k     ^  ^  -■ 


(-W  ^30-^^1^ 


g(Xg,y) 


[1-0 (i)]  . 


The  asymptotic  behavior  for  h(k)   is  determined  similarly; 

(1^9)  h(k) ^  [1   +  0(i) 

Since  w(y)   increases  monotonically  we  have 
w  -  w  =  ^^;(Y  )  -  \-j{y)    >  0  for  all  y,  0  5  Y  <  Y  .  and  hence, 
from  ([|.8),  g(y,k)   increases  exponentially  with  k  for 
Re  k  >  0,   But  we  knox-7  from  (25)  that   g  dies  aviay  like  t- 
as   k  — >  +CO  ♦   It  follows  that  each  coefficient  in  T„(Y  ,k) 
must  vanish.   It  can  be  proved,  see  [l^]  ,    that  the  successive 
vanishing  of  the  coefficients  in  T„(Y  ,k)   yields  a  series 
of  equations  relating  the  derivatives   P(y)   and  Q(y)   at 
y  =  y   and  these  relations  through  (1+)  and  (5),  determine 
G  (X  ,y)   in  terms  of  G(X  ,y). 


•iv  G/'    m-;    '?;    ;.  • 


(X.^}0   ^ 


xxo    cvij^rq«'(j:;3a   r-' r    ,0.1    - 


ca-:.v    -in.j"vi': 


-!     'v  \-l-v ;  >jL 


\  .""  /    •'     —  t 


[y/fd     'yc: 


/  ,f 


'    .1      ni    ;• 


i^    ir:c^;. 


•~'     '  V  i ', 


.,  1 1  ■  •■    "    ■• 


\  ■  ■/  ■  I 


oi 


^\\V\S>Z 


o  ■  -  "i; 
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For  the  present  we  simply  state  that   P(y)   is  such  that 
(50)  e^T^(Y^,k)  -0 


for  Re  k  >  0. 

From  equation  (k-Q)   we  then  have 

G(X^,y)    / 


(51) 


^(^o*y^  (        \\  .  ,     - 

g(y,k)-  1 [l   +   0(i)j  for    |k|    >  kp 


Rek>k,    0<y<Y, 
p*         -    J    -      o 

For  Y   <  y  <  0  we  find  exactly  the  same  behavior, 
oo  ~ 


Part  III.   Asymptotic  Behavior  of  the  Flow  Near  the  Center* 

11.   The  Inverse  of   g(y,k). 

We  have  assumed  that  the  two  integrals, 

(^   -k(x-X^) 

g(y,k)  =  j    e        G(x,y)dx 

X 
o 

r^  -k(x-x  ) 

h(k)  =  J    e      °  H(x)dx 
^o 

converge  absolutely  for  Re  k  >  f'^  .   V/e  assume  the  solution 

G(x,y),  H(x)   is  continuous  diff erentiable  on  X  <  x  <  oo , 

0  <  y  <  Y  .   3y  the  general  theory  for  the  inverse  of  a 

Laplace  transform,  we  have  then: 

r  %'^^'^ 
(1)  lim        2?fl  J  e^^^^~-S^   g(y,k)dk  =  G(x,y)        x   >  X^ 

np — >co  K  -iT 

o 

,^  +iT      ,   ,      V    ^ 
o  k(x-X    ) 

11m       ^\  e  °  h(k)dk  =  H(x)  •     X   >  X^. 


'i'    nav: 


(V 


.'^   %^ 


V>    (,-.-    '  V 


ov^r:'    a^nct    3'><'    \ 


.0 


cr 


:-[   <    !>J 


<-    I 


;>::yja    ynJ 


oo 


,  :.iit,v.)^ 


-  V  .bli 


(    X-x 


cr.), 


y         1    ^-    ( 


CO 


■ch{yj  U 


[   -  w;;:-i 


, 'T^  >  X  >    ^X     no   sic:;;;..  Jao^j11J:j:    u:,;> 


f  - 


<    X  1, 


38. 

By  shifting  the  path  of  integration  in  these  integrals 

from  Re  (k)  =  K^      to   Re  k  =  r'       where  K-,    <     ;<'  ,  it  is 

o  11    "-o 

possible  to  determine  the  asymptotic  behavior  of   G(x,y)  and 
H(x)   as   X  — >  00  • 

We  first  investigate  the  continuation  of  g(y,k)   and 
h(k)   for  Re  k  <  K^.   From  (20-2i;),  we  see  that  g(y,k)   and 
h(k)   are  defined  in  the  left  half  k-plane  except  at  (a)  the 
singularities  of  g-j_(y,k),  g-(y,k)   and  g.  (y,k)   and  (b)  the 

zeros  of   D(k)  =  g^Si^j  -   "^i^^j   +  (Pl+'^q]_  )gl. 

¥e   assume  for  the  present  that  the  singularities  of 
g(y,k)   and  h(k)   are  all  simple  poles  at  the  points   k  =  k 
with  residues   g(y,k  )   and  h(k  )   respectively.   ^Je  shall 
study  them  in  greater  detail  later. 

Prom  Part  II,  equations  (1+9)  and  {$!) ,    we  have   |kh(k)| 

and   |kg(y,k)|   bounded  for   |k|  >  k   Re  k  >  k   where   k 

may  be  chosen  arbitrarily.   By  a  well-kno^-m  application  of  the 

Tauberian  theorem,  the  behavior  of   G(x,y)   as  x  — — >  oo   is 

given  by 

. —  k  (x-X  ) 
(2)  G(x,y)  r.)        e  ^    °   g(y,k^,) 


n 

and  that  of   H(x)   by 

•r k  (x-X  ) 

(3)  H(x)  -}__  e  ^    °  hCk^). 

n 

Whether  G(x,y)   and  H(x)   die  out  as  x  — *>  co   depends  on 


^'>' 


has,    '         '         V:.   -toiv.v;-;o!jf 
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whether  any  of  the  quantities   k   have  positive  real  parts. 

We  now  investigate  the  singularities  of  g(y,k)   in  order 
to  determine  these  quantities   k  ,   Singularities  in  g(y,k) 
may  occur  if  D(k)  =0  or  if  g^(y,k),  g^(y,k),  g|^(y,k)   are 
singular.   We  will  show  first  that  the  singularities  of  g,, 
g-,  g,   do  not  affect   g(y,k)   and  h(k).   Then  we  will  examine 
the  values  of  k  for  which  D(k)  =  0, 

IS,  The  Singularities  of  g-|_(y,k),  g^(y,k),  g,  (y,k). 

As  was  pointed  out  in  Section  (9),  the  functions   g, (y,k), 

g-(y,k)   and  g,  (y,k)   are  defined  for  all  k  ^Inless 

a^  +  b  k  =  q  where   q  =  0,  -1,  -2,  •••,  etc.,  to  k  =  a 
00  -i.  J         3         1  1  J  q 

q-a  ■ 

=  — r —  •   a   is  always  negative  since   a   >  0,  b   >  0.   It 

o 
can  be  shown  now  that  g,  ,  g-,   and  gi   all  have  simple  poles 

at   k  =  a„   and  that  the  residues  are  all  proportional  to  the 

singular  solution  g2(y*k).   Thus  for  k  — >  a   where 

a„  +  b  a  =  q,  we  have 
o    o  q   ^* 

(14.)   g^(y,k)(k-a^)  --  Jiqg2(y*^q^  +  (k-aq)<{)3_(y,k) 

g^(y,k)(k-a^)  -v.  j^^g2(y,a^)  +  (k-a^)(|)^  (y,k) 

g^(y,k)(K-a^)  ^  j^qg2(y,aq)  +  (k-aq)(l)^(y,k) 

where  (}),(y,k),  <|)^(y,k),  ij),  (y,k)  are  appropriate  functions 
analytic  in  the  neighborhood  of  k  =  a  and  j,  ,  j'  and 
jj    are  appropriate  constants. 
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On  substituting  these  expressions  in  T(y,k)   and  using 
(ij.)  once  more,  we  find  that   (k-a  )T(y,k)   is  analytic  near 
k  =  a  .   Thus  both  T(y,k)   and  D(k)   have  simple  poles  at 
k  =  a  ,   Accordingly,  from  (5),  g(y,k)   is  not  singular  at 
k  =  a  .   Similarly  we  find  that  h(k)   is  not  singular  at 

M. 

k  =  a  . 

q 

This  means  that  the  functions   g(y,k)   and  h(k)   can  have 
singularities  in  the  k-plane  only  for  those  values  of  k  for 

which  the  denominator  D(k)  =  g^S^.Y   "  ^1^3y  "*"  (Pi'*"^^l^S]_ 
vanishes.   Note  that  this  expression  i^  independent  of  the 
particular  choice  of  g^» 


13.   The  Eigenvalues. 

¥e  proceed  next  to  examine  the  roots  of  D(k) 
=  I^iij  -   SiS3y  +  (Pi+kq^)gl«   f'or  large  values  of   |ki   and 
Re  k  >  k   we  have,  see  ikh)    Part  II, 

D(k)  ~  k  f^io^'^'^ao  ^  "ii^  "^  "^^Iv 

so  that  D(k)  has  no  roots  in  any  right  half  plane  if  |k| 
is  sufficiently  large.  In  other  l■^fords,  D(k)  has  at  most  a 
finite  number  of  roots  with  Re  k  >  0, 

Since  any  two  regular  solutions  of  dL(g)  =  >/ +  mk  +  nk 
differ  only  by  a  multiple  of  g^(y,k),  it  follows  that 

(9)     D(k)  =  i^g^y  -  g3ygi  +  (Pi+kqi)ii 
=  (g'tiy  -  gyii)  +  (Pi+kqi)gi 

=  (-gjj  +  Pi  +  kq^)ii 
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where   g""(y,k)   is  the  regular  solution  of  ^(g")  =  j^+  mk  +  nk 

for  which  g'"(Y^,k)  =  g"""  =  0.  Thus   D(k)  =0  if  Fy  =  P^  +  ^^2.* 

In  other  xirords,   D(k)  =  0  for  the  eigenvalues   k   of 
the  following  problem: 

(10)  ^(g'"  )  =  /  +  mk  +  nk^ 

g"(Y^,k)  =  0 
gy(Y^,k)  =  P3_  +  kq^ 
g"(y,k)   regular  at  y  =  0. 

In  order  to  determine  the  residue  "gij,^    )      of  g(y,k) 
at   k  =  k  ,  we  determine   g(y,k)   in  terms  of  g"(y,k)   instead 
of   g,(y,k).   From  (20)  -  (2[|.)  of  Part  II,  we  find  then, 

(11)  g(y,k)  =  .  Jt  g  +  1  ^y  _ Ji  ly o       ^    ^  ^    ^ 

gl  Sl^-gy  +  Pi  +  kq^] 

We  assume  that  ~  (g^)  f^  q-^  for  k  =  k^..  It  follows  then 
from  (11)  that  the  residue  g(y,k^)  of  g(y,k)  at  k^  is 
given  by: 

(12)  g(y,k  )  =    ^^l  .    i  %, —  \  s  ij,^^) 

Sii"  ak  Sy  •*•  ^1/ 

=  r  g""'(y,k  ). 

Similarly  the  residue   h(k^)   of  h(k)   at   k  =  k^  is  foiond 
to  be 

(13)  h(y,k^)  -  n^. 

Two  eigenvalues  of  k  are  easily  determined,  namely 
k  =  0  and  k  =  X,   We  prove  this  first  in  the  case   k  =  0. 
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From  the  equations  for  the  basic  flow,  (lli)  Part  I,  we 
see  that  if   U^(n),  C^(:n),  P^Cn)   is  a  solution  then  ^^i^^.U+tJ ) , 
C^iMl+t^)) ,    ^oillC^-^^Q^)      is  also  a  solution.   Expanding  U^, 
C^,  P^   in  powers  of   e^  we  find  by  (llj.),  Part  I, 


P^(^(l+e^))  =?^{rj)   -.  e^P^ 


^0(7)  (-^o  >) 


This  corresponds  to  the  special  solution  of  (2l|.  and  (39), 
Part  I,  that  is  independent  of   r: 

A 


(15) 


Pl(r|,r)  =  J^irp   =^ 
Ci(>],r)  =  C,(n)  =  ^ 


'1^^ 


l^r  =5" 


H^(r)    =  ^H  =  -1  . 


Corresponding  to  this  solution,  equations  (5B)  and  (60) 
Part  I  for  G(x,y)  and  H(x)   have  the  solution 


(16)  G(x,y)  =  ^G(y)  =  -  ^ 
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H(x)  =  qH  =  -1. 


Thus    G(y)   satisfies,  see  (58),  the  equation 
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that  is    G(y)   is  a  solution  of  (10)  with  k  =  0. 

That   k  =  X   is  also  an  eigenvalue  follows  from  the  fact 
that  if   u(r,t),  c(r,t),  p(r,t)   is  a  solution  of  the  equations 
for  spherical  flow  (1)  Part  I,  then  u(r,t+e-|_),  c(r,t+e-|^), 
p(r,t+e^)   is  also  a  solution.   This  means  in  the  present  case 


^"-^^     "  =  t^  Il^(r-Nt.e^)),  c  .  ^   0^(r-Nt+e^)), 


P   = 


2  2 
°'  ^   2   ^o^^"'^^*'^^!^^   ^^  ^  solution  of  (1)  Part  I  where 


-X 


(t+e^) 


U  (ri),  C  (n),  P  (n)   is  Guderley's  fundamental  solution. 
Expanding  in  powers  of   e-,   according  to  the  procedures  of 
section  i|  we  find  that 
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C^(r,r|)  =  x^l^^'7^  " 
P^(r,>|)  =  xPi(r,>|)  = 
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is  a  solution  of  (2l+)  and  (39)  Part  I  if  we  take  H^(r)  =  -r~  /H^. 
This  solution  corresponds  to 


G  =  ^G(x,y) 


X(x-y) 


X-1 


X(x-Y^) 
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H  =  ^H(x)  =  - 


H. 


being  a  solution  of  (08)  and  (60)  Part  I  f jr  G(x,y)   and 
H(x).   Substituting  (I8)  in  (08)  Part  I  we  find  that 
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(19)  .g(y)  =  H  ® 
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satisfies  X(,g)  =  ^  +  mX  +  nX^,  ,g(Y  )  =  0  and  ,g  (Y  ) 
/V  K         o  A.  y   o 

=  P]_  +  ^q]_»   Thus   T,g(y)   is  an  elgenf unction  and  X   an 
eigenvalue  of  the  system  (10), 

The  other  eigenvalues  are  not  readily  calculated.   In 
appendix  II  a  rather  crude  method  for  calculating  k   and 
g" (y,k  )   is  described.   It  appears  from  these  preliminary 
calculations  that  the  only  non-negative  values  of  k  are   0 
and  X  and  that  the  highest  negative  value  of  k  is  approx- 
imately -2» 

It  can  also  be  shown  that  D(k)   vanishes  at  least  once 
q-a  q+l-a 

between   a  =  —r —  and  0.   ^-,    =   — r — -     if  q   is  a  sufficiently 

o         ^  o 

large  negative  Integer  so  that  an  infinite  number  of  eigenvalues 

exist. 


II4..   Behavior  of  the  Flow  Quantities  at  the  Center, 

Substituting  equations  (I3)  and  (II4.)  in  the  asymptotic 

expressions  (2)  and  (3)  for   G(x,y)   and  H(x),  we  obtain 

the  asymptotic  representations: 

r —  k  (x-X  ) 
(20)         G(x,y)  -  )   e  ^    °  T  ^^  s"' (j ,1^^) 
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for  0  <  y  <  Y  ,   The  sarae  behaviour  can  also  be  shown  for 
-  ''  -   o 

Y   <  y  <  0, 
oo  -  "^  - 

Using  (lt|)  in  the  expressions  for  V,,  ¥„,  and  V-,  (equa- 
tions (53),  (55),  (56)  and  (i+O)  Part  I)  we  obtain 

r k  (x-X^) 

(21)  V^(x,y)-\_e^    ""    V ^     Sy^^^^n^ 

n 

V2(x,y)  -^  e""^  """"^   r^  [p2+V2-(Vl^2)s''(y,k^)] 

i, ji 

n 

V3(x,y)  --^^  e  ^  ""  °  r^  [p3-(V^^3^s'''(y,k^)] 
n 

or  using  (27)  Part  I  and  (29)  Part  I, 

r ~k 

(22)  U^(r|,r)-)_r  T^     '^{^\*K'^ 


n        ^ 


r -k 


Pl(<|,r)-2_^  '"k^Pl^^^^n^ 
n 

Ci(n,r)  -.y~  r  ^T^  <^i"(^^k  ) 

'      ^ n     ' 


where 


k  (y-X  )   „ 
(23)      ^i(^»^n^  =  ^  ""    °  Sy^^'^n^ 

Pl(V],k^)  =  e^"  ^"  °   (P2+V2-^V^2^s''(y'^n^^ 

k  (y-X  )  ,. 

c'^'(y),k^)  =  e  ^   °  (P3-(k^-tt3)g"  (y,k^)). 

Since   k  =  0  and  k  =  X   are  two  of  the  values  of  k^ 
it  would  appear  from  (22)  that  the  flow  quantities   U^,  P^ 
and  C,   do  not  die  out  as   r  — >  0,   However,  if  we  examine 
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the  original  formulation  of  the  problem  in  sections  I],   and  5> 

Part  I,  we  see  that  we  are  always  free  to  choose  for  any  actual 

physical  situation  the  txrro  constants   R,   and  t",   (see  (3l(-) 

Part  I)  which  fix  H(X  )   through  (51)  and  H^(X  )   through 

(52),   We  recall  that   R   +  eR,   is  the  position  from  which 

the  shock  starts  and  t   +  et,   is  the  time  at  which  it 

o     1 

starts. 

We  could  set  about  proving  that  it  is  possible  to  choose 
R,   and  t,   in  such  a  way  that   P    and  P  ^      in  (22)  vanish 
identically  and  hence  that   U,  ,  C   and   P.,   all  vanish  as 
r  —t-O. 

However,  we  proceed  to  eliminate  these  roots  in  a  more 
specific  fashion.   Prom  section  II4.,  pages  53  and  Sk,    we  have 
-  =  t%  U^(r-Nt.s,)),  p  =  ^  P^(r-Nt.e,)), 


c  =  :^  C^(r"''(t+e^)).   But  u  =  ^^  U^(r  t(l+e^)), 

p  =  ££  P  (r"'*'t(l+e  )),  c  =  ^  C  (r"'^t(l+e  ))   are  two  solutions 
^    o  O  U   o  o 

of  the  equations  of  flow  (l)  Part  I.   Here   e^  and  t-^     are 

considered  of  the  same  order  as   e,  see  section  ij..   Combining 

them  we  see  that 


(2li) 


is  a  solution  to  (1)  Part  I,   Expanding  in  powers  of   e^  and 
e^  and  using  (15-18),  (10)  and  (23)  we  obtain 
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(25)   X  I  u(r,t)    =  U^(r|)    +  e^u'"(rj,0)    +   e^r"^  u''(k|,X) 

+  lower  order  terms   in      e    ,    e, 

o'   1 

X  ^  c(r,t)  =  C^(vj)[l  +  e^c"""(y,0)  +  e^r'-^c^Ct^  ,X) 

+  lower  order  terms  in   e  ,  e,  ] 

X  ^  p(r,t)  =  PQ(n)[l  +  eoP'"'(n,0)  +  e^r'^p''''(^^,X) 

+  lower  order  terms  in  e  ,  e,  ] • 

Now  from  (22)  and  equation  (23)  Part  I  we  have 

-k 


(26)  U(r|,r)  =  U^(t^)  +  eU^(r,vi)  --U^(>j)  +  e^  ^  Tj,  """^^'^n^ 

-k 
C()],r)/C^.,()^)  =  1  +  eC^(r,r|)->l  +  e^  r  ""Pi,  c"(V],k^) 

-k 
P(r|,r)/P^(t^)  =  1  +  eP3^(r,r|)'^  "^  ""  ^T"  ^  ""^k  P'"(^»kn^* 

If  we   choose      e      =   e  P'       and     e,    =   e  P..  ,   we   obtain,   using 
o  o  1  A. 

(2i|.-26)   and     n  =  r~^t, 

(27)  U(r|,r)'-^  ^;:Y-^  U^  ((l+e^)  (n+e^r"^))    +  eyr     ""r  ^  ^'''(l»l^) 

'  -'•''"^1*1     r  \  /  ^ —  n 

-k 


^(^'^^^    ,,      ll   -X   ^oO^-^^o^r^l^'^O    "   'P"  ""^k  <^"(^'^n^ 
l+e^n     r  ^  /  ^ —  n 

l+e-,r)     r  ^  '  ' —  n 


Here   k   no  longer  assumes  the  values   0  and  X. 

It  is  clear  then  that   U(Y|,r),  C()^,r)   and  P(n,r)   differ 
less  and  less  from  one  of  Guderley's  solutions  as   r  — >  0 
since  the  values  of   k   are  all  negative. 
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Appendix,   Computation  of  the  Elg:envalues» 

In  section  I3  we  required  the  eigenvalues  of  the  following 
non-homogeneous  system:   g(y,k)   is  a  function  which  satisfies 
the  non-homogeneous  equation: 

(1)  ^(g)  =  yg"  +  (a+bk)g'  +  (c+dk+ek^)g  =  i+mk+nk^ 

see  (16)  Part  III,   The  coefficients   a,  b,  c,  d,  e ,  Jl ,   m,  n 
are  all  analytic  functions  of  y. 

The  initial  conditions  satisfied  by  g(y,k)   are 

(2)  g'(YQ,k)  =  Pi  +  kq^ 

(3)  g(YQ,k)  =  0. 

In  addition,  for  special  values  of   k,  g(y,k)   is  regular 
at  y  =  0,   These  special  values  of   k  are  the  eigenvalues  we 
seek.   We  would  like  to  approximate  them  n'umerlcally  and  show 
that  there  are  no  eigenvalues  of  k  with  positive  real  part 
except   k  =  A  =  1.39^7* 

We  reduce  the  problem  to  one  which  Is  easier  to  handle 
numerically.   Suppose  that  k  is  one  of  tlae  eigenvalues. 
Setting  g(y,k)  =  g^(y,k)il;(y,k)   where   g^(y,k)   is  the  regular 
solution  of  (^(g-,  )  =  0  and  g^  (0,k)  =  1,  see  (11)  and  (lii.) 
Part  II,  we  find  using  (1-3),  the  differential  equation  for 

^(y,k) 

ih)  ygl^"  +  (2yg{  +  (a+bk)gi)\l/'  =  j2,  +  mk  +  nk^ 

and  the  initial  conditions 
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k9, 

(5)  vl/(Y^,k)  =0  at  y  =  Y^ 

(6)  gj_il/  +  g^il/'  =  P-L  +  kq^     at   y  =  Y^ 
or  using  (5) 

(7)  gjL^'  =?!•*■  ^^1* 

Prom  the  condition  that  g(y,k)  and  g, (y,k)  are  regular 
at  y  =  0 ,  \<!e  see  that  il/(y,k)  and  hence  \1/'  (y,k)  are  regular 
at  y  =  0,   Setting 


2  ^o^^^'     r   ^-%^(b-b^)k 


(8)  ^(y,k)  =  g^^y  "  "  expj   ^ ^:^  dy 

0 


we  find  using  {I4.)   and  (7), 


(9)   \l/(y,k)vl/'(y,k)  =  J   y^^  vKy^^klgJ^^ly  ,k)(i(y^) 

Y.  ■"• 


+  m(y3_)k  +  n(y^)k^)dy^  +  (p^+kq3_)g]^^(Y^,k)4/(Y^,k) 


or 

y 


(10)  i|/'(y,k}  =  ^~^(y,k)     y-^vl/(y^,k)g~^(y^,k)U(y^) 

^0 

+  m(y^)k  +  n(y^)k^)dy3_  +  C(k)F-'-(y,k) 


where  y 

~  o 


r'^T.f^r         V^o."l/•-<r    K'\(P(tt      \+rr^(^T      \V+-n  ( ^T  Wc^  ' 


(11)   C(k)  =-J   y^^v|<(y-L,k)g--^(y^,k)(i(y^)+m(y^)k+n(y3_)k^)dy^ 
0 

+  (p-L+kq^)g"^(Y^,k)?(Y^,k). 
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Using  (8)  in  the  expression  (10)  for  \|/' (y,k)  we  see  that 
\V'(y,k)  is  regular  at  y  =  0  only  if  C(k)  vanishes,  that  is 
from  (11), 

(12)  (vj_+^<i2.^el-^i^^,'^)^r^o'^^   =     y'^^(y,k)g^^(y,k) 

UiY-^}    +  m{j^)k  +   n(y^)k2)dy^. 

Eqiiation  (12)  is  therefore  a  condition  on  the  eigenvalues 
of  k. 

Setting 

T  a^+b^k 

(13)  gl  (y,k)\l/(7,k)  =  y  °  °  F(y,k) 

we  rewrite  (12)  as 

Y 
a  +b  k  f  °  a  +b  k-l 

(llj.)   (p^+kq^)Y^°   °  F(Y^,k)  =  J    y  °   °   P(y,k) 

0 

(J?(y)+in(y)k+n(y)k2)dy. 

From  (8),  (13)  and  ^(g-|_)  =  0  we  obtain  the  differential 
Bquation  for  F(y,k),  namely 

(15)         yF"  +  (a'''"+b"'''k)F'  +  (c'"'+d""'k+e'""k^)F  =  0 


where 


a'"'(y)  =  -a(y)  +  2a(0) 
b"(y)  =  -b(y)  +  2b (0) 

c'"'(y)  =  -(a(y)-a(0))y"^(a(0)-l)  -  a'(y)  +  c(y) 
d'''(y)  =  -(a(y)-a(0))y-^b(0)  -  (b(y )-b(0)  )y"^(a(0)-l) 
-b'(y)  +  d(y) 
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e'''(y)  =  -(b(y)-b(0))y~^b(0)  +  e(y). 
From  (8)  and  (I3) 

(16)  P(y,k)   is  regular  at  y  =  0 
and,  using  g, (0,k)  =  1 

(17)  P(0,k)  =  1. 

The  method  for  calculating  the  values  of  k  is  now 
as  follows.  We  expand  P(y,k)  in  increasing  powers  of     y 
using  (15-17): 

(18)  P(y,k)  =  1  +  P'^(k)y  +  P2(k)y^  +  ...  . 

We  also  expand  ^(y)  +  m(y)k  +  n(y)k   in  powers  of  y: 

(19)  /(y)  +  m(y)k  +  n(y)k2  =  I^  +  \W7   +  l2(k)y^  +  •••  • 

Substituting  (I8)  and  (19)  in  (llj.)  and  integrating  we 

obtain, 

CD         00     yn     n_ 

(20)  (p, +kq.  )  ]         P  y"  =  ^    — TvA-r-  Y   F  I 

n=0        n=0  °   °    m=0 

To  find  the  first  approximation  for  the  first  two  eigen- 
values we  take  just  the  first  terms  in  the  sums  in  (20). 
Prom  P  =  P(0,k)  =  1   and   I^  = /(O)  +m(0)k  +  n(0)k^  x^e 


obtain 

"o  ~o"' 


(21)       P   +  kq   =  Ag)^^(0)ktn(0)k^ 
V^-LV       P-L  -^  i^q^        a  +b  k 


or 
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(22)  P^(k)  =  (p-,_+kq^)(aQ+b^k)  -  (/(0)+ra(0)k+n(0)k^)  =  0. 
On  substituting  the  appropriate  numerical  values  we  find 

(23)  p-'-(k)  =  ,1776  +  .1^683k  +  3.3960k^  =  0, 
Solving  (23)  for  k  we  obtain 

k  =  -.037     or     1.1^16* 

For  the  second  approximation  we  use,  see  {!$), 

(0)  \      P  -   c'"'(0)+d'''(0)k+e'""rnlk^  =  .  c'"'(0)+d''''(0)k+e'"'(0)k^ 
^        a"(0)+b"(0)k  %'^''o^ 

along  with 

(25)  I-^  =  P  (0)  +  m'  (0)k  +  n'  (O)k^, 

The  relation  corresponding  to  (20)  can  then  be  reduced 
to  the  polynomial 

(26)  P^(k)  =  (aQ+bQk)P^(k)  +  y^[ (a' (0)+b' (0)k) (a^+b^k) 

-  c^  -  d^k  -  e^k2][P^(k)  -  yi(P3_+kq^)]  +  y3_(a^+b^k) 
W  (0)  +  m»(0)k  +  n'  (O)k^)  =  0 

or  substituting  numerical  values, 

P^(k)  =  .018250  +  5. 10815k  -  l.[^25l39k^ 

+  l,623692k^  -  .191l86k^  . 

This  method  can  be  continued  indefinitely.   The  approx- 
imation to  (20)  yields  successively  various  rational  relations 
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53. 

to  be  solved  for  k.   Each  rational  relation  reduces  to  a  poly- 
nomial whose  degree  Increases  by  two  at  each  step. 

Thus,  as  it  stands,  this  method  appears  to  yield  two  new 
eigenvalues  at  every  step.   The  computations  were  made  however 
by  rejecting  the  higher  order  terms  in  k  and  dealing  with 

polynomials  which  increase  by  only  one  degree  in  k  at  each 

2 
step.   Thus  we  reject  the  fourth  order  term  in  k  in  P  ,  the 

fifth  and  higher  terms  in  P-^   and  so  on.  Each  approximation 

then  yields  one  nev;  eigenvalue. 

The  third  approximation  comes  then  from  the  poljmomial 

p3(k)  =  +.1226  +  10.2li.90k  -  ,31369k^ 

-  l4..2308k^  -  ,52936k^ 

The  approximate  eigenvalues  are  given  in  Table  I,  The 
eigenvalues  appear  to  be  all  negative  except  k  =  0  and 
k  =  A  =  1,395«   The  first  negative  eigenvalue  is  approximately 
-2. 


Table   I, 

^1 

^o 

^-1 

l.i|l6 

-.037 

1.390 

-.0087 

-2c3 

l.i^lO 

-.012 

-1.8 

lo395 

0 

„ 

^-2 

P^(k) 

P^(k) 

P^(k)  I.I^IO     -.012     -1.8       -7.6 

Correct  value 
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